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We provide analytic expressions for the Green's functions in position-frequency space as well as for 
the tunneling density of states of various Luttinger liquids at zero temperature: the standard spinless 
and spinful Luttinger liquids, the helical Luttinger liquid at the edge of a topological insulator, and 
the Luttinger liquid that appears either together with an ordering transition of nuclear spins in a 
one-dimensional conductor, or in spin-orbit split quantum wires in an external magnetic field. The 
latter system is often used to mimic a helical Luttinger liquid, yet we show here that it exhibits 
significantly different response functions and, to discriminate, we call it the spiral Luttinger liquid. 
We give fully analytic results for the tunneling density of state of all the Luttinger liquids as well 
as for most of the Green's functions. The remaining Green's functions are expressed by simple 
convolution integrals between analytic results. 

PACS numbers: 71.10.Pm, 75.30.-m,71.70.Ej, 85.75.-d 
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I. INTRODUCTION 

Helical conductors have received much attention re- 
cently. They consist of one-dimensional (ID) conduc- 
tors exhibiting spin-filtered transport, for instance, of 
right- moving modes carrying spin up, (i? f), and left- 
moving modes carrying spin down, (L ],). Such con- 
ductors appear on the edges of two-dimensional (2D) 
topological insulators (where spin should more correctly 
be addressed as Kramers partner) in nanowires made 
of tree-dimensional (3D) topological insulators^— in ID 
semiconductor wires with strong spin-orbit interaction 
in an external magnetic field (with approximate spin- 
filtering) ,— ~— or in carbon nanotubes, where the combi- 
nation of spin-orbit interaction and strong external elec- 
tric fields can cause helical conduction] 14 ' 15 Perfect spin- 
filtered conduction modes can also appear in regular ID 
conductors in the presence of nuclear spins as explained 
in further detail below. Aside from the spin-filtering, 
if a helical conductor is brought in the proximity of a 
superconductor it supports Majorana bound states at its 
endsf 1 ^— Such Majorana end states have attracted much 
interest recently as parts of elementary excitations with 
nonabelian statistics that may be useful for topological 
quantum comput at ion . 2 3 ' 2 5 ~ — 

While much of the discussion of these systems in the 
literature is focused on the fundamental physics and the 
possible applications of the helical conductors, it is still a 
challenge to provide a unique experimental proof of the 
existence of the helical states. In this paper we show that 
probing the spectral properties of ID conductors provides 
an extra tool to investigate experimentally the differences 
between helical and non-helical systems. This comple- 
ments recent proposals to probe spin-dependent trans- 



port properties in helical conductors using, for instance, 
a quantum-point-contact setup or a Scanning Tunneling 
Microscope (STM).22r^ 

Using the Luttinger liquid (LL) theory, which allows us 
to take into account the effects of the electron-electron in- 
teractions, we compute the position-frequency dependent 
Green's functions for the helical conductors proposed in 
the literature and, for comparison, for the regular spin- 
less and spinful LLs. This allows us to obtain informa- 
tion about the spectral properties of these systems. In 
particular, we give explicit expressions for the local (or 
tunneling) density of states (DOS) that may be probed 
directly, for instance, by STM experiments. Most of our 
results are fully analytic which, to our knowledge, has 
not been achieved before, even for regular LLs. We show 
that electron-electron interactions play a central role and 
determine much of the spectral properties. In particular, 
we show that the interactions cause a substantially dif- 
ferent response of the conduction modes in ID quantum 
wires as compared with the conduction modes at the edge 
of a topological insulator. 

Furthermore, our results provide a handle to detect nu- 
clear spin order in a ID conductor. As mentioned above, 
it has recently been shown that the presence of nuclear 
spins in an electron liquid can lead to new physics at 
low-enough temperatures. In a semiconductor, such as 
GaAs, nuclear spins and electrons interact via the hy- 
pcrfinc interaction. This generates a long-ranged inter- 
action among the nuclear spins mediated by the elec- 
tron gas, which is of Ruderman-Kittel-Kasuya-Yosida 
(RKKY) type. In 2D and ID conductors, electron- 
electron interactions increase the long-range nature of the 
RKKY interaction (for the ID case, see Ref. |34| ) and can 
trigger a finite temperature ordering of the nuclear spins 
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in systems of finite size<2£r— The nuclear magnetic order- 
ing has only little effect on the electron state in 2D ) 35 ' 36 
which may hinder a possible detection of this order via 
electronic transport or spectroscopy. 

The situation differs in fD metallic systems, where 
electron-electron interactions have a dramatic effect. The 
Fermi liquid description breaks down and the physics 
is dominated by collective low energy excitations. The 
appropriate description of this scenario is known as 
LL theory (see, for instance, Refs. and |4l| for re- 
views). Signatures of the LL physics have been ob- 
served nowadays in many systems such as in semicon- 
ducting GaAs quantum wires ^r— in metallic single wall 
carbon nanotubes^ - — in fractional quantum Hall edge 
states^ - in bundles of NbSc3 nanowires^i in polymer 
nanofibers^ in MoSe nanowires^ 3 - in conjugated poly- 
mers at high carrier densities,— and in atomically con- 
trolled chains of gold atoms on Ge surfaces^ 

On a spinoff from the LL theory, it has been shown in 
Refs. [37] and HH that in the presence of a lattice of nu- 
clear moments and of the hypcrfine interaction, an exotic 
ordered phase emerges at low temperature in which the 
magnetic moments of the nuclei and the conduction elec- 
trons of a ID conductor are bound together due to the 
RKKY interaction. The nuclear spins order in form of a 
helimagnet. In turn, the resulting magnetic Overhauser 
field strongly affects the conduction electrons, much in 
contrast to the 2D case. It induces a gap for one half 
of the low-energy modes, which yields a partial electron 
spin polarization that follows the nuclear spiral polariza- 
tion. The rest of the low-energy modes remains gaplcss 
and can be described by an effective LL model with per- 
fectly spin-filtered left- and right-moving modes. We call 
this LL state a "spiral Luttinger liquid" (SLL). A simi- 
lar SLL model has been shown to describe the quantum 
wires with spin-orbit interaction in the presence of an ex- 
ternal magnetic field, in which the spin-dependent band 
shift induced by the spin-orbit interaction is commensu- 
rate with the Fermi momentum kpr^ The latter type of 
LL corresponds thus to a special point in the general LL 
description of interacting quantum wires with spin-orbit 
interaction^ - — 

The SLL should not be confused with the helical Lut- 
tinger liquid (HLL), which has been introduced to de- 
scribe the edge states of 2D topological insulators^ The 
HLL is also spin-filtered, but differs from the SLL, the 
most important difference being that the SLL describes 
the bulk properties of a ID system, while the HLL de- 
scribes the edge of a 2D system, or the surface of a 
nanowirc made of a 3D topological insulator^ For non- 
interacting electrons, the low-energy physics of HLLs 
and SLLs is equivalent and, as mentioned, this equiv- 
alence has been of much interest very recently because 
both HLLs and SLLs support the Majorana bound states 
at their ends if they are brought in the proximity of a 
superconductor! 1 ^— For interacting electrons, the prox- 
imity induced superconductivity and the shape of the 
Majorana end states is subject to strong rcnormaliza- 



tion and can deviate strongly from the noninteracting 
Hmiti^Ir— Here we add a further clement of caution: 
We show that the assumed equivalence of HLL and SLL 
breaks down for any interacting electron system, and that 
major differences arise in the response functions. These 
differences persist if the systems are brought in the prox- 
imity to a superconductor, and may lead to a different 
response of the bulk as well as of the Majorana properties 
in both types of systems. 

Because of the specific magnetic ordering in the SLL, 
we find that its electronic excitations can be divided into 
a regular sector that exhibits a typical LL (or HLL) type 
behavior, as well as an irregular sector that exhibits a 
more special type of physics. In the non-interacting limit 
the quasiparticles in the regular sector are typical chiral- 
propagating modes, while the irregular sector is gapped. 
There is still some propagation in this sector because 
the left-moving spin-up (L f) an d right-moving spin- 
down {R 4-) modes, characterized by the boson fields 4>l^ 
and 4>Rli are n °t separately gapped; instead, they are 
pinned together through (j)^ = 4>Ri- ln this sense, the 
gapped sector exhibits similarities to the LL edge states, 
for which the reflection of the quasiparticles at the end 
of the wire pins the L and R movers together i£r£L The 
difference between the two is that the pinning of 4>l-\ an d 
4>fli for the irregular sector of the SLL is a bulk prop- 
erty, while for the edge of the LL the pinning of the L 
and R modes is local and happens only at the end of the 
wire. Thus, the modification of power-law exponents in 
the SLL irregular sector is similar to the modification of 
the end-tunneling exponents for a LL. In contrast, this 
physics is absent in a HLL. 

In this paper, we analyze the spectral properties of a 
SLL in comparison with those of the regular LLs and 
the HLL. In particular, we analyze the dependence of 
the two-point retarded Green's function on position and 
frequency, as well as the dependence of the DOS on fre- 
quency 

Our first result is a closed and simple analytical 
form for the frequency dependence of the local retarded 
Green's functions for the spinless LL, the HLL, as well as 
for the regular sector of the SLL. While the analytical ex- 
pressions of the Green's function of a LL in position-time 
[x, t) space as well as momentum-frequency (k, uS) space 
are well known ; 40 ' 41 ' 82 " — the position-frequency [x, lo) 
dependence has been either calculated directly numeri- 
cally by Fourier transform 8 ^ or has been obtained by an 
efficient recursion method^ However, to our knowledge, 
an explicit analytic expression has not been provided yet. 
Here we succeed in directly evaluating the Fourier trans- 
form of the position-time dependent Green's function, 
and we obtain a closed compact form for the position- 
frequency dependence of the imaginary time and retarded 
Green's functions. For the spinful LL as well as for the ir- 
regular sector of the SLL we cannot provide such a closed 
form. Instead, we give the results in form of simple con- 
volution integrals. 

Second, we use the obtained results to derive the fre- 
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quency dependent local (or tunneling) DOS. We compare 
the results for the SLL to that of the standard LL and 
of the HLL. For HLLs and regular LLs the irregular con- 
tributions are absent. The DOS does not show a gap 
and follows the usual power-law dependence. The DOS 
of the SLL is given by the sum of the contributions for 
the regular sector and the gapped irregular sector and 
therefore exhibits a pseudogap. Inside the pseudogap 
the DOS shows a power-law dependence, characteristic 
of the regular sector. However, at the gap edge, the 
DOS is dominated by the irregular sector and depends 
on the value of the interaction parameter, such that for 
weak interactions it exhibits a typical divergence, while 
for strong interactions this divergence goes to zero. Con- 
sequently, we propose to use the dependence of the DOS 
on frequency to test the apparition of a magnetically or- 
dered phase due either to the hyperfine interaction with 
the lattice of nuclear spins or to the spin-orbit interac- 
tions. The measurement of the DOS can provide a direct 
access to the size of the pseudogap, while the form of the 
gap-edge singularity may allow to infer the strength of 
the interactions. A similar behavior of the DOS can be 
found in ID systems with charge density wave order £ — 

We note, finally, that our work has been followcd-up 
by Ref. IqoL which contains a computation of the spectral 
properties of the SLL in (fe, w) space, and a derivation of 
the DOS near the pseudogap edge beyond the harmonic 
approximation used in the present work. 

The plan of the paper is as follows: In Sec. HH we 
introduce the various models under consideration, cor- 
responding to the standard LL, the HLL, and the SLL. 
Section IIIII constitutes the core of the paper. In this 
section we present our results for the two-point (x, uj) 
Green's functions. In Sec. IIV1 we compute the DOS of 
the various LLs. We provide a summary of our results in 
Sec. El Three appendices are also included, containing 
details of the Green's functions calculations (Appendix 
fA"]) . the DOS calculations (Appendix [B]), as well as the 
derivation of an important identity (Appendix [Cjl . 



II. MODELS 

We start by introducing the models that we will con- 
sider in our analysis. These models consist of different 
variations of the LL model. We summarize here only the 
main features of these models, and for a thorough discus- 
sion of the LL physics we refer the reader, for example, 
to Refs. and [ij]. Our notations and normalizations 
for the bosonic fields follow the conventions of Ref. |4l|. 
We consider the following three types of LLs: the stan- 
dard spinless and spinful LL; the helical LL appearing on 
the edge of topological insulators; and the spiral LL ob- 
tained by the self-ordering transition due to the hyperfine 
coupling to nuclear spins, or due to a commcnsurability 
transition in systems with spin-orbit interactions. 



A. Standard Luttinger liquids 

Interacting spinless electrons in a ID conductor can be 
described by the LL Hamiltonian 



H 



LL 



dx 



K 



(y<j>) 2 + K{vef 



(i) 



where </> and 9 are bosonic fields such that —Vip/n mea- 
sures the electron density fluctuations in the system and 
V0/7T is canonically conjugate to (f>. We set h = 1 
throughout the paper. The x integral is performed over 
the system length C. We assume C to be much longer 
than the width of the low-energy wave packets, CkF <C 1, 
where kp is the Fermi momentum, which allows us to use 
a continuum description. The LL parameter K incor- 
porates the effects of the electron-electron interactions, 
K = 1 corresponds to non-interacting electrons, and 
< K < 1 corresponds to repulsive electronic inter- 
actions. The quantity v is the velocity of the bosonic 
model, which for an ideal LL is given by v = vf/K, with 
vf being the Fermi velocity. 

In this theory, the electron operator ip is written as 
"0 = i-'L + ipR, where ipL corresponds to the left-moving 
modes with momenta close to — kp, and ipR to the right- 
moving modes with momenta close to kp . These fermion 
operators are related to the boson fields as 



ip r (x) 



Vr 



(2) 



with r — L, R = — , 



and with the chiral boson fields 



4> r {x) = 4>{x) — r9(x). 



(3) 



The quantity a in Eq. ([2]) is a short distance cutoff, 
limited from below by the lattice constant, and r\ r is the 
ladder operator, or Klein factor, whose effect is to lower 
the r particle number by 1. 

For electrons with spin, the Hamiltonian of the ID 
conductor in the LL regime takes the form H = H c + H s , 
where H c s have the same form as Hll in Eq. (JTJ) , 



dx 
— 'i 
2tt 



(V<^) 2 +K v (V6„f 



(4) 



for v = c, s labeling the charge and spin degrees of free- 
dom, and v v = vf/K u . The decoupling of H c and H s 
is the consequence of the spin-charge separation in ID. 
Here K c s = 1 corresponds to the nonintcracting case, 

< K c < 1 to repulsive electron interactions, and K s ^ 1 
to a broken spin SU(2) symmetry. The presented LL the- 
ory applies only for K s > 1 , as a value of K s smaller than 

1 would lead to the opening of a spin-gap and so to dif- 
ferent physics. The normalizations are chosen such that 
— V0 C;S v / 2/t'" give the magnitude of the charge- and spin- 
density fluctuations. The electron operators can again be 
separated into left- and right-movers, ipo- = ipLa + i^Ra, 
where 



lp r <r( x ) 



nra_ irkpx —ir4> rl7 (x) 



(5) 
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with the additional spin index a — J \,\ r — +, — and the 
boson fields 

4>ra{x) = -j=[(p c {x) - r9 c {x) + a(4> s {x) ~ r9 s {x))]. (6) 

This expansion of <p ra (x) contains both charge and spin 
fields, indicating that an electron is a superposition of 
both spin and charge degrees of freedom. 



B. Helical Luttinger liquid 

The quantum spin-Hall effect appears in various sys- 
tems with strong spin-orbit interaction and preserved 
time-reversal symmetry. While the band structure has 
been long known^^Ir— such systems have found much re- 
newed interest recently because it was recognized that 
the electron state indeed describes the quantum spin- 
Hall effect, a topological state of matter characterized by 
a bulk gap and gapless edgesi^ Recent experiments on 
HgTc quantum wells have provided direct evidence for 
nonlocal transport in the quantum spin-Hall regime in 
the absence of any external magnetic field, in agreement 
with the theory 94 ! 95 Due to the strong spin-orbit cou- 
pling, the left- and the right-movers of the gapless edge 
modes carry opposite spin. These edge chiral excitations 
are described by a separate class of one-dimensional LLs, 
the HLLs i 70 ' 96 Contrary to the chiral LLs at the edges of 
fractional quantum Hall systems, these edges states do 
not break time-reversal symmetry. 

The linearized Hamiltonian of the non-interacting HLL 
in the fermionic language can be written as 



H HLL 



-Vf 



dx {ily m id x tp m - ipl^dxipLi). (7) 



One can easily describe the interactions between fermions 
of the same species, as well as the forward scattering 
processes in the bosonized language i 30-32 ' 70 ' 97 ' 98 Using 
Eq. ([3]) the Hamiltonian in Eq. ([7]) can be written in a 
bosonized form. Moreover, in a HLL, the spin and chiral- 
ity indices coincide due to the strong spin-orbit coupling, 
and, as one can see directly in Eq. ([7]), the spin index 
is redundant. It is therefore more convenient to treat 
the Hamiltonian of the HLL as an effectively spinless 
LL by introducing the fields $ = (<pu^ + 4>l\)/^ and 
G = {4>li — </\Rf)/\/2. When including the interactions, 
the bosonized form of the Hamiltonian of the interacting 
HLL becomes 3 ^- 3 ^?^ 



1 t V 



1 



K 



HLL 



-(v$) 2 + /f ffLL (ve) 5 



(8) 



where we have introduced the parameter Khll, which is 
a generalized LL parameter taking into account the in- 
teractions, and v — vf/Khll, the velocity of the collec- 
tive excitations. Recent Quantum Monte Carlo simula- 
tions have shown that this HLL description of the helical 



conductor remains indeed valid as long as the interac- 
tion strength remains smaller than the band gap pro- 
tecting the edge states from the higher-dimensional bulk 
states^^ 



C. Spiral Luttinger liquid 

The presence of nuclear spins in a ID conductor can 
substantially modify the LL behavior. It was shown in 
Rcfs. [33 and HH that the hyperfine interaction between 
the nuclear and electron spins can indeed trigger a strong 
feedback reaction between nuclear spins and the electron 
modes, leading to ordered phases in both subsystems. 
This order is manifested in the apparition of a nuclear 
helimagnet that is stabilized by its coupling to the elec- 
tron system. In turn, the resulting spiral nuclear mag- 
netic field generates a relevant sine-Gordon perturbation 
for the electrons. The Hamiltonian that captures this 
electron physics can be written a o 37 i 38 



B eS 



■ cos 



with 



v+ = (v c K c + v s K s l )/K 
v- = (vcK- 1 +v s K c )/K, 



where we have introduced 



K = k c + k: 



(9) 



(10) 
(11) 



(12) 



In Eq. ([9]) we have neglected unimportant marginal 
terms coupling the "+" and "— " fields. 

The field B e g in Eq. (JfJ]) represents the effective Over- 
hauser magnetic field induced by the nuclear spiral order. 
The fields </>± and 6± are related to the charge and spin 
boson fields by^ 8 .^ 



(13) 




(14) 

(15) 
(16) 



An identical sine-Gordon renormalization was also iden- 
tified in low-density ID electron conductors with spin- 
orbit interaction, such as GaAs or InAs wires,— and can 
also appear in carbon nanotubes . 14 ' 15 Promising can- 
didates are furthermore Gc/Si core/shell nanowires as 
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they combine an unusually large, controllable spin-orbit 
interaction ; 102 ' 103 a good hole carrier confinement in the 
core, and high mobilities ,i2i~ifll 

The Overhauser field is then replaced in the semicon- 
ductor wires by a uniform magnetic field, and in the car- 
bon nanotubes by the combination of intrinsic spin-orbit 
splitting together with an electric field. The relevant 
rcnormalization occurs as a commensurability transition 
when the electron density is tuned such that kp = k so , 
with k so being the inverse spin-orbit length. 

The cosine term in Eq. ^ is strongly relevant and the 
field 0+ becomes pinned to a constant. Excitations in 
this sector are characterized by a gap A. The nuclear 
Overhauser field binds indeed one half of the electron 
modes into a spin and charge mixing electron density 
wave with a spiral spin polarization that follows the nu- 
clear spin helix. Through the hyperfine interaction the 
stability of the nuclear helix is further enhanced by the 
SLL. The combined clcctron-nuclear-spin order was pre- 
dicted to be stable up to temperatures of 10-100 mK 
for experimentally available GaAs quantum wires^ r 44 ' 108 
and 13 C substituted single- wall carbon nanotubesj 1 ^-!^ 

The remaining half of the electron modes can be de- 
scribed by a modified LL, characterized by a different 
LL parameter. We have termed the resulting low-energy 
Hamiltonian the "spiral Luttinger Liquid" (SLL). Simi- 
larly to the HLL discussed in Sec. Ill Bt the SLL possesses 
spin-filtered left- and right-movers, yet the SLL exhibits 
important differences with respect to the HLL. We will 
discuss these differences in detail in the following. 

Let us mention that Hamiltonian (|9|) does not directly 
apply to the case carbon nanotubes, as the latter require 
a two-band description that captures the physics of the 
two Dirac valleys of the nanotubesj 113 ' 114 However, as 
was shown in Ref. HH, cosine terms such as in Eq. © 
open gaps in each Dirac valley separately and overrule 
the otherwise strong inter- valley coupling caused by the 
electron-electron interaction. The SLL physics of carbon 
nanotubes can be, therefore, captured by two indepen- 
dent copies of Eq. ©, one for each valley. In the fol- 
lowing, therefore, we restrict entirely to the single-band 
model of Eq. ([9]) and note that for most of the measurable 
response functions the difference to nanotubes amounts 
in a mere factor 2 in the results. 



D. Differences between Spiral and Helical 
Luttinger liquids 

The main difference between the SLL and the HLL is 
the origin of the spin-filtered conduction modes. While 
for the HLL they are a consequence of the special 2D 
or 3D band structure, for the SLL they result from the 
opening of a spin-dependent pseudogap A in an oth- 
erwise conventional spinful ID bulk conductor. In the 
non-interacting limit the separation between the gapless, 
spin- filtered modes (L I), (R f) and the gapped modes 
(L f), (R I) is perfect, and for energies |w| <C A the 



physics of the SLL, governed by (L J.), (R f), is identical 
to the HLL. For interacting electrons, however, parts of 
the gapless phase, described by the boson modes </>_ and 
9- affect the (L f)j (R D sector as well and lead to what 
we call the "irregular" response even at low energies. 

To see this in detail, let us project the initial electron 
Hamiltonian onto the gapless ("— ") sector and rewrite 
the Hamiltonian §§§ in the standard form of a spinlcss 
LL |T|), reminiscent of the HLL. We start with rescaling 
the fields d>_ , 6L as 



2K ° A 

~kk~ h 



KK 

-2K7 6 -' 



(17) 



which allows us to rewrite the gapless sector of Eq. 
as 



dx 

H SLL — I -^VSLL 



1 



■<y<j>¥ + K SLL (vof 



(18) 



_KsLL 

with the parameters 

K SLL = 2K C /(K C K S + 1), v SLL =v_. (19) 



If we neglect the "+" fields for the time being, the boson 
fields are projected onto the gapless sector as 



KK S ' 

Vs K 1 

and the left- and right-movers become 



<!>li 

0Lt 



A K C 



(20) 
(21) 
(22) 



(23) 
(24) 

(25) 



where K = K c - Kj 1 . 

While the form of 4>l\. and <pg^ is similar to that in 
Eq. (|3]), the presence of the propagating 8 field in the 
expression of (j>m and 4>l\ is, at first glance, surprising, 
since the exponentials of these boson fields correspond 
to the nonpropagating fermion operators of the pinned 
spiral density wave. However, such a single-electron pic- 
ture is inappropriate for an interacting system. Instead, 
for an interacting system, a left- moving electron cor- 
responds to a superposition of left-moving particlc-like 
and right-moving hole-like bosonic modes that together 
carry unit charge and spin. Only in the noninteracting 
limit does the amplitude of the right-moving excitations 
vanish. Moreover, the gap A builds up from the elemen- 
tary bosonic modes, which carry only a fraction of the 
electron charge. The connection between cf)^ and (j)^ 
is therefore a natural consequence of electronic interac- 
tions. A similar connection holds between <fin^ and 4>li- 
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Since <j>ni = — </>Lt (after projection), only the antisym- 
metric part 9 between the gapless (j)^ and <pL^ can be 
excited, which is encoded in Eq. (|25|) . 

This peculiar 9 dependence is the major difference be- 
tween the SLL and the HLL, and leads to the irregular 
contributions to the response functions that are detailed 
below. It is a pure interaction effect and, indeed, the 
prefactor K of the irregular 9 correlators of the fields in 
Eq. (|25[) vanishes in the non-interacting limit. 



A. Standard Luttinger liquids 

1. Spinless Luttinger liquid 

For the standard LL, the calculation of imaginary-time 
fermionic correlators reduces to calculating the correla- 
tors between the boson fields <f> and 9. To this end, it 
is convenient to introduce left- and right-moving eigen- 
modes 



III. GREEN'S FUNCTIONS IN (x,oj) SPACE 



rKt 



(32) 



This section describes the derivation of one of the main 
results of the paper, the form of the two-point Green's 
functions in (x, w) space, which constitutes the basis 
for the description of the spectral single-particle prop- 
erties of the conductors. While much is known about 
the Green's functions in position-time (x, t) (and imag- 
inary time t = it) space ; 40 ' 41 as well as in momentum- 
frequency (fc,w) space ; 82 ' 83 to our knowledge the corre- 
lators in the (x,w) space have been to date determined 
only numerically . 85 ' 86 Here we provide analytical, closed 
expressions for the spinless LL, the HLL, and the regular 
sector of the SLL Green's functions. For the spinful LL 
as well as for the irregular sector of the SLL we provide 
integral formulas in form of convolution integrals. We 
focus here on imaginary-time ordered Green's functions 
and on real-time retarded Green's functions, yet our cal- 
culations can be easily repeated for all other types of 
Green's function. 

We define in the standard way the imaginary time 
(Matsubara) Green's function in (x, t = it) space for 
a translationally invariant stationary system as 



G T ,f'{x,T) = -(T T V4(Z,T)VV(0,0)>, 



(26) 



where f = (r, er, . . . ) describes the chirality, spin, and any 
other quantum number characterizing the electron oper- 
ators, and T T is the usual time order operator. The real 
time Green's functions are expressed using the elemen- 
tary greater and lesser functions 



G^ fl (x, t) 



-i{4{x,t)iPr>(0,0)), (27) 
-i(Vv-'(0,0)v4(M)>, (28) 



which allow us, for instance, to express the retarded 
Green's function as 



Gl%{x,t) = #(t)[Gl f ,(x,t) - Gf tF ,(x,t)], 



(29) 



with iD(t) the unit step function. The Fourier transforms 
to frequency space are given by 



G" r,(x, uj) 



— oo 

oo 



dre^" T G p ,f(x,t), (30) 
dte luJt G^,{x,t) 7 (31) 



where a e {>,<,ret}. We restrict our analysis to zero 
temperature. 



which commute with each other. The evaluation of 
their correlators in (a;, r) space is well described in 
textbooks ! 40 ' 41 At zero temperature, for a translation- 
ally invariant system described by the Hamiltonian in 
Eq. ((J), we have 



([Mx,T)-4> r (0,0)} 2 ) = - In 



x — irvT 



ira 



(33) 



This result is valid in the limit |ar|,w|r| 3> a, in which 
bosonization applies. The electron correlators are then 
evaluated using the relation 



/ e i[\<t>(x,T)+»8(x,T)] e -i[\<t,(0,0)+»8(0,0)]\ 



e 2 



j([A0(z,t)+/x0(:j;,t)-A0(O,O)-/j0(O,O)] 2 ) 



(34) 



for arbitrary constants A and /i. Using Eqs. (|2|) and (|3|), 
for the spinless LL this leads to 



G r ,r(x, 7") 



2ira 



ir r 



Ml 



1 It 



X + ivr 



with 



7r = 7-r = 



K + K- 1 - 2r 



(35) 



(36) 



Note that — 7r| = 1- Cross-correlators G r . r > with 
r =/= r' vanish as they do not preserve the number of lcft- 
and right-movers. 

The Fourier transformation to G rir (x, iuj n ) can dir ectly 
be evaluated using standard integral tables [Ref. Ill5l 
3.384.9], and leads to a Whittaker function. Yet since the 
properties of the Whittaker function might be unfamil- 
iar to the reader, we propose here a modified calculation 
that expresses the result in terms of the more commonly 
known Bessel functions, in a combination that is equiva- 
lent to the Whittaker function. 

Since this calculation yields one of the main results of 
this paper, we provide here the details for the imaginary 
time Green's function calculation. The analogous calcu- 
lation for the retarded Green's function can be found in 
Appendix IaI As explained below, the correspondence be- 
tween these two Green's functions, usually given by the 
Wick rotation, is actually problematic here (and would 
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be even more obscure if we would use the Whittaker func- 
tion) , and an independent calculation for the two Green's 
functions is necessary. 

We start by rewriting the imaginary-time Green's func- 
tion as 



' *p <p ^X 5 



47T7 



■ {d rx H~ 3iVT ) 



rx — ivt rx + ivt 

(37) 

with 7 = j R = (K + K'- 1 - 2)/4. We use differenti- 
ation to reduce the exponent of the singularity in the 
integrand, which provides us with an elegant way to reg- 
ularize the unphysical ultraviolet divergences. Indeed, we 
have < 7 < 1 for K > 0.17, which will assure conver- 
gence of the integrals below. For smaller values of K, we 
can use an analytic continuation of the result in 7 and 
neglect any occurring divergence, or we can use again the 
rcgularization by further applications of the differentia- 
tion. Wc thus need to calculate the integral 



I(x,iu) n ) 



x 2 + {vt)' 2 



(38) 



whi ch can be done using standard integral tables [Rcf. 
Til 8.432.5]. The result is given by 



I(x,i(J n ) 



2ay/ir ( 2\x\v 



2-7 



JCidzwnlM (39) 



where K a (z) is the modified Bessel function and T(z) 
Euler's Gamma function. To obtain the Green's function 
we need to evaluate also the derivatives in Eq. (|37]l. The 
time derivative yields a factor —uj n /v, while the deriva- 
tive with respect to rx is evaluated using the identity 



— z- a K Q {z) = z~ a K a+l {z), 
az 



(40) 



from which we obtain the result 



G ri r(x, ioj n ) 



_ e irk F x iund ( 2\x\v 



2^r( 7 + l) v 2 \\u n \a 2 
K y _i(\xu n \/v) - sign(rx(J n )K y+ i (\xu) n \/v) 



(41) 



This expression of the Green's function holds for a real 
uj n . The analytic continuation to the retarded Green's 
function by the Wick rotation is, however, not well de- 
fined, due to the nonanalytic dependence on |w n |, which 
arises because taking the zero temperature limit and the 
analytic continuation do not commute. We circumvent 
this problem by directly computing G ret (x 1 iS) using con- 
tour integrations. The details are summarized in Ap- 
pendix [AJ The result is 



G r r et r {x,uj) 



Lj + a f 2i\x\v 



20Fr( 7 + l) v 2 



K r _i(|x|w + /i«)-sign(ra)iC Y+ i(|x|a; + /iu) , (42) 



with uj + = Lu+iO. As a function of w, the Bcsscl functions 
have branch cuts along the negative imaginary axis, i.e., 
at the Matsubara frequencies uj n < 0. This shows that 
indeed the analytic continuation between the real and 
imaginary time Green's functions via the Wick rotation 
cannot be naively used, justifying the independent calcu- 
lation. 



2. Luttmger liquid with spin 

The LL model for true, spinful electrons can be ob- 
tained by considering two copies of the LL model (for 
the charge and spin fields), as described by the Hamil- 
tonians in Eq. (j4]). The electron operators are given by 
Eq. (0, together with Eq. (j6]). Since the charge boson 
fields commute with the spin boson fields, the evaluation 
of the bosonic correlators as expressed in Eq. (f34| factor- 
izes into two exponentials, each of them depending only 
on one of the charge and spin fields. Thus Eq. (|3"5"|) can 
be replaced by 



Gra,ra 



2ira 



1 7»v 



X — IV C T 



HI 



x + iv c r 



—ia 


It, 8 


ia 




x — iv s T 




x + iv s T 





with 



1r,v T— r,u 



K v + K- 1 - 2r 



, (43) 



(44) 



where v = c,s. As before, the cross-correlators with 
(r, a) 7^ {r',a') vanish, and \"f r ,v — l'rv\ = 1/2- Generally 
we have v c 7^ v s , and the Fourier transform cannot be 
calculated using Eq. (f3"5)l . To our knowledge, a general 
solution to this Fourier transform has not yet been found. 
However, it can be evaluated exactly in the special cases 
of v c = v s and of x = 0. If the velocities coincide, the 
Green's functions can be obtained by the same type of 
calculation as the Green's functions for a spinless LL, and 
are given by Eqs. (|4Tj) and (|42)l , on the replacements r — > 
(r,a) and 7 -> j r , c + j r , s = [(K c + K s ) + {R- 1 + K^ 1 ) ~ 
4r]/8. The limit x = 0, important for the calculation of 
the DOS in Scc. lIVl can be obtained in a similar manner, 
noting that one needs to evaluate an integral of the same 
type as the one in Eq. Q38p , on the additional replacement 
v — > ^/v c v s , and on taking the limit x — > 0. 

For general velocities and for x 7^ 0, we note that Eq. 
(|4"3"]) can be viewed as the product of two spinless Green's 
functions with labels c and s. Hence its Fourier transform 
can be written as the convolution 



G 



r (x,iuj n ) 



(2\x\/a) 



l-7c-7s 



4nv c v s T( 7c + l)r( 7s + 1) 

X / d^'nfr,c{ X ' iu n)frA X ' iu; n (45) 
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with 



fr tV (x,lU) n ) = 



<->v 



Und I \oj n \a 
V 



V 



I' 



x K lij _i{\xu> n \/v u ) - sign(rxu n )K 7v+ i (| xoj n \/v„) 

(46) 

While Eqs. (|4"5]l and (j4*u]) do not provide direct insight 
into the analytic shape of the Green's functions of the 
spinful LL, due to the exponential convergence of the 
integrals they are convenient for numerical evaluation. 
A similar expression holds for G^* rcr (x, w), yet with a 
slower convergence because of the power-law decay of the 
Bessel function along the imaginary axis. 



B. Helical Luttinger liquid 

The helical edge states arise from the special band 
structure of a topological insulator. They are described 
by two counterpropagating edge modes with opposite 
spin (Kramers doublets), say (L |) and (R f)- In the 
non-interacting system, an electron with spin j" tunnel- 
ing into the system has to move to the right and can- 
not enter into a coherent superposition of (R f) and the 
(nonexisting) (L f). In the interacting system, however, 
chiral electrons are no longer the good eigenstates and 
the tunneling electron excites a coherent superposition 
of particle-like excitations of (L f) and hole-like excita- 
tions of (R I), together carrying unit charge and unit spin 
f. The resulting theory is a fully regular LL as described 
by the Hamiltonian in Eq. (|5|). The Green's functions 
are consequently given by Eqs. (|4"Tj) and (|4"2j) with the 



replacement 7 ->• jhll = (Khll + K 



HLL 



2)/4. 



C. Spiral Luttinger liquid 



The SLL corresponds to the low-energy theory of con- 
duction electrons in the presence of the pseudogap A. 
Contrary to the HLL, the number of degrees of freedom 
is doubled, with one gapless sector and one gapped sec- 
tor. If the SLL is caused by the nuclear spin ordering, 
there are two intrinsic energy scales in the system, the 
nuclear spin ordering temperature T* and the pseudogap 
A of the electrons, fulfilling for most of the semiconduc- 
tor quantum wires ksT* < A, with ks the Boltzmann 
constant. We will assume henceforth that the tempera- 
ture T satisfies T < T* such that the SLL state remains 
stable. We will also consider only frequencies \ui\ ^> fc^T, 
at which our zero-temperature analysis remains valid. 
For the SLL in spin-orbit interaction split semiconduc- 
tor quantum wires, only the condition \u\ ^> ksT is 
imposed as the scale T* is absent. The distribution of 
energy scales allows us thus in particular to explore the 
regime \u\ ~ A, which will be of much interest below. 

In order to calculate the Green's function G ra - tra -(x,T) 
in imaginary time, we first use Eqs. (O and ^ to write 



the electron operators in terms of the charge and spin 
bosonic fields, and then use the change of basis given 
by the set of Eqs. (|T5|) — (fT5 |) to write them in terms of 
(f>- 1 9^ 1 <j) +1 + . After these manipulations, the Green's 
function takes the form 



G /■ 



2?ra 



-F-(x,t)F+(x,t), (47) 



where 




-(l + <TrW-=2(M*,T)-M0 J 0)) 

V s 



+ (-S_+a^K^\ (M*,t)-M0,0)) 



(48) 



and 




lK\ (*e-g)(Ms,T)-^ + (0,0)) 



(a-r)(e + (x,r)-e+(0,0)) 



(49) 



1. Regular sector 

Let us first consider the Green functions associated 
with the modes (L 1) and (R f), which corresponds to 
ra = 1. In this case, only the correlator ((4> + (x,t) — 
(f>+(0, 0)) 2 ) enters in . Since the field </>+ is gapped, 
((<f) + (x, t) — </>+(0, 0)) 2 ) decays exponentially on the scale 
\J x 2 + v 2 t 2 ~ v/A, which implies that F 1 + (x,r) w 1. 
The correlation functions between the modes (L !) and 
(R f) are therefore described by a standard spinless LL 
theory. This corresponds to the "regular sector" of the 
SLL. The Green's functions of the regular sector are given 
by Eqs. ([4"T]) and (|4^|). on the replacements r = L — > (L l 
), r = R — > (R t), and the new parameters v — > vsll — 
V- as defined in Eq. pT|) and 



7 -> 1; 



L SLL 



K SLL 



reg 



(50) 



with Ksll = 2K C /{K C K S + 1). Cross-correlators be- 
tween the (L l) and (R f) modes vanish. 



2. Irregular sector 

The correlation functions between the modes (L \) or 
(R I) correspond to ra = — 1. The expression for these 
correlators is more complicated since F^ 1 in Eq. (|47[) 
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can no longer be approximated by unity due to the pres- 
ence of the 8 + field in Eq. (|^|) . which is conjugated to 
the pinned </>+ field and so strongly fluctuates. We de- 
note this as the "irregular sector" since the correlation 
functions cannot be inferred from a standard spinless LL 
analysis. 

The Green's function Gf f(x,iui n ), with f £ 
{(Lf), (iJJ,)}, requires the calculation of the Fourier 
transform of the product FZ 1 (x,t)F^ 1 (x,t), for which 
we unfortunately we have not been able to find an analyt- 
ical, closed form. We can, however, evaluate F^-y(x,iu} n ) 

individually and express Gr,r{x,iu) n ) in form of again a 
convolution integral similar to Eq. (|45|) . In addition, 
we will see in the next section that the knowledge of 
F^(x = 0,iu) n ) is sufficient to obtain a closed form for 
the DOS. Therefore, we will focus here on the individual 
evaluations of F^(x, ioJ n ). While an analytic expression 
can be found for FZ 1 (x, iu) n ) for all x, this is possible for 
F^ ± (x, iuj n ) only at x = (as shown in the next section). 
For i^Owc provide the result in form of a convolution 
similar to Eq. (|4"5)) . 

Let us start with F^ 1 (x 7 iuj n ). In order to compute the 
asymptotic behavior of the associated Green functions, 
one nee ds to make a number of approximations. Follow- 
ing Ref. Ill6l we first replace the sine-Gordon term in Eq. 
© by its quadratic expansion i.e. B e g cos(V2A'0 + ) — s- 
BeffK^. The calculation o f F^ becomes then similar 
to th e one per form ed in Ref. Ill6l . The result is (see also 
Refs. [Il7| and [nit) 



F+ (i,t) we 



2 ^2 



(51) 

where m 2 = B c gK/av+, K = K c + 1/K S as defined be- 
fore, and K+ is the strong coupling LL constant of the 
"+" sect or of the model defined by Eq. ([9]) . Contrary to 
Ref. Ill6l K*_ is only little renormalized by the renormal- 
ization group flow of the cosine term in Eq. (j9]),— and 



we can use K T 



1. 



Since the quadratic expansion of the cosine term in 
Eq. ([§]) may seem a rather crude approximation, one 
can check that Eq. (p)Tj) makes sense by comparing it to 
the non-interacting result obtained for K c = K s = 1. 
In this limit the cosine term in Eq. (|9]) is equal to 
cos(20 + ), and corresponds to the backscattering of R i 
and L f fermions (the ones associated with the "+" sec- 
tor). Moreover F^{ m (x, r) is proportional to the Green's 
function of these massive fermions and therefore can be 
written as 



duj n dq tuJnT+lqx -iu n -vq 
) 2 w 2 + A 2 + v 2 q 2 ' 

(52) 

with A cx yj B c gv I a and v = v + = v_ . At x — 0, we can 
evaluate F^[ m (0,T) analytically and find 



F+r(0,r) 



A 

2^ 



*i(A|r|), 



(53) 



where Ki is the modified Bcsscl function. Since K i(z) ~ 
c~ z /y/z for z> 1, the asymptotic behavior of Fti(0,r) 
in the non-interacting limit is consistent with the behav- 
ior described in Eq. (|5ip . This indicates that the rather 
crude approximation used to calculate F^ 1 (x,t) for a 
given K captures correctly the asymptotic limit for the 
Green's function. 

The full (x, iui n ) dependence of FZ\ cannot be de- 
termined analytically. We can, however, give a so- 
lution in terms of a convolution integral, similar to 
Eq. (gSJ), if we write Ft x (x,r) = fi(x, r)/ m (x, r), 
with fi the irregular part from the gapless sector, 
fi(x, t) = [a 2 /(x 2 + v 2 _t 2 )] 1 / 2K , and f m the massive part, 

f m {x, t) = e~^y /x2+v l T2 . The gapless part fo(x, iu„) is 
identical to I(x,iuj n ) in Eq. (j35]) with the replacements 
7 — > 1/2 A' and v — > v + and the solution given by Eq. 
(1391) . The t integral for the massive part f m {x, iuj n ) can 



be evaluated through standard integral tables (Ref. Ill5l 
3.914) and we find 



f m (x,iu) n ) = 2 / dr cos(u} n r) e" 
Jo 

" K x (\x\P/v+), 



Zk\f x2 + V 



(54) 



with (3 = 



4K 2 



u} 2 . Using the latter expressions, 



F+ 1 (x,iu„)= daj' n fi(x,iuj' n )f m (x,icj n -ioj' n ), (55) 



which is again well suited for numerical evaluation be- 
cause of the exponential convergence of the Bcssel func- 
tion K\ . 

The second factor for the Green's function in Eq. (|47l) , 
FZi(x, t), is entirely evaluated within the gapless LL sec- 
tor and is given by 



FZ^r) 
with 



-sign(r) 



—ia 


firr- 


ia 




x — iv-T 




x + iV-T 





K 2 



li 



K 2 K S K S {K C — K7 1 ) 2 



AK SLL K 2 8K C K 8K C (K C + AT 1 ) ' 



(56) 



(57) 



The sign(r) factor in Eq. (|5rJj) is necessary to ensure the 
correct fermion statistics, which is no longer enforced by 
the form of the power-law factor, since it depends now 
only on x 2 + vtr 2 (the sign can also be verified by an 
analytic continuation of the independently calculated real 
time Green's function). In contrast to the regular LL, 
both exponents in Eq. (|56j) arc identical, expressing the 
fact that the left- and right-movers are pinned to each 
other [see Eq. (|2"5j) ]. 

The Fourier transform resembles closely the integral 
I(x,iuj n ) of Eq. (|38]l. with the supplementary factor 
sign(r) in the integrand. The result can a gain be eval- 
uated using standard integral tables [Ref. Ill5l 3.771.1, 
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Ref. Ill9l . 12.2.3]. It is less compact than for I(x,iuj n ) 
and given by 



f_i(i,«w n ) = -isign(w„) 



7 7irr-i(i a;a, ™l/' i; -) 



a^/nT(l - 7 irr ) 
2v- 



2uj n a 2 



XV- 



y irr (\xUn\/V-) 



(58) 



where L a (z) is the modified Struve function and I a (z) 
the modified Bessel function. 

We again perform an independent calculation for the 
retarded Green's function associated with FZi, denoted 
by FZ{ ret , for which the sign sign(r) — > sign(t) leads to 
the further terms described at the end of Appendix [A] 
The result is 



for the cross-correlators is identical to the result of the 
auto-correlators determined above, on the replacement 
x —> x~ x' everywhere, except in the phase factor e lrfcira: , 
where we make the change x — > x + x' . 



IV. DENSITY OF STATES 

The DOS per unit volume a (or tunneling DOS) p(w), 
can be obtained from the imaginary part of the retarded 
local Green's function as 



P (oj) = --lmG ret (x ->■ 0,lj), 

7T 



(60) 



where G ret = £ G r r f r<T . 



F_{ [x, u) = -sign(o; 



1 



V- I r(7 4rr ) \ Lu + a 



2i\x\v- 



x Ki irr Qx\u + /iv-) 



T(l - 7irr) 



2xv- 



H i _^ i (\u J x\/v-)-Y i _(\uJx\/v-) 



(59) 



where H a (x) denotes the Struve function, Y a (x) the 
Bessel function of the second kind, and u + = u) + iO. 
The part depending on the Bessel function Ki_^. rr cor- 
responds to I(x,uj) given in Eq. (|A5[) . while the other 
terms are parts of the integral that cancel for I(x, w) but 
contribute here because of the additional factor sign(i). 
They are given by Eq. (|A7|) . We notice, however, that 
the latter terms are all purely real so mainly the com- 
plex Ki _ 7irr dependent term contributes to the spectral 
properties of the system. The pro perties of the Struve 
functions are described in Ref. Ill9l . Sec. 12. The Struve 
functions are related to the modified Struve functions 
by ~L a (z) = — ic _1 2" Q H Q (iz). The asymptotic behav- 
ior of the Green's functions at large z is determined 
by noting that Hi_ 7 (z) — Yi_ ry (z) ~ 0(z~2~ 7 ) and 

I 7 _i(z) - Li_ 7 (z) - 0(>-3-t), so Eq. dSS]), as well 
as the last terms in Eq. (|59|) have the large-distance 
non-oscillating decays | | 1 and |x|~ 1 |w|~ 27irr , re- 
spectively. For z — > 0, the Struve functions go to zero, 
while -L,_i and Yi_„ have divergences that arc regular- 

'221 

ized by the ultraviolet cutoff. 

The convolution of the F^(x, iu> n ) times the prefactor 
— e~ lrkpX /2ira yields the Green's function G rr7 . rr7 (x, r). 

Let us, finally, note that cross-correlators f ^= f' 
in the irregular sector, with f , f ' e {(L f), (R |)}, 
are nonzero. The corresponding Green's functions 
(T T ipf(x,T)ipl,(x',T')) actually depend on the two posi- 
tion variables x and x'. This is because, for the cross- 
correlators, the response is no longer translationally in- 
variant due to the phase factors e ±lk r( x + x ). This spa- 
tial variation is the only difference between the auto- 
correlators and the cross-correlators. Hence the result 



A. Density of states for the Luttinger liquid and 
the Helical Luttinger liquid 

The calculation of the DOS for a spinless LL is detailed 
in Appendix [51 Starting from Eq. (|42|) one obtains 



p(uj) 



2 1 



'27rur(l + 27)' 



(61) 



with 7 = (K + K^ 1 — 2)/4 and A a = v/a on the order 
of the bandwidth. The prefactor 2 arises from the sum 
over r ~ L, R. Here we have dropped a divergent, fre- 
quency independent contribution to the DOS. The latter 
is a well-known artifact from the LL theory in the short 
distance limit x — > 0, and must in practice be limited by 
the value of the short-distance cutoff a. The form of the 
DOS in Eq. (foTj) corresponds to the results known in the 
literature ! 40 ' 41 For non-interacting electrons, 7—^0 and 
the DOS becomes a constant, reflecting the fact that a 
linear dispersion is assumed in LL theory. 

On the replacement K — > Khll in the expression of 
7, the result (fBTj) holds as well for the HLL. 

For a spinful LL, we have noticed in Sec. IIII A 21 that 
the x — >• limit can be evaluated using the Green's func- 
tions of a spinless LL while replacing 7 — > ~/ c + j s = 
(K c + K~ x + K s + K- 1 - 4)/8 and v ^/vJh. Since 
v cs = vf/K C:S , we can incorporate this rescaling of veloc- 
ities by redefining A a = ^/v c v s /a = vp j a\J K C K S . This 
yields 



p{u) = 4 



| w/Aa |2(7e+7 s 



27™r(l + 2(7 c + 7s ))' 



(62) 



where the four channels (r, er) contribute identically to 
the DOS. 



B. Density of states for a Spiral Luttinger liquid 

For a SLL, the DOS is the sum of the contributions 
from the regular and irregular sectors. The spin-resolved 



11 



DOS is 



A» = —Im[G r £(x -> 0,w) + G^(z -> 0,w)]. (63) 

Independently of er, one of the two contributions in 
Eq. (jB"3"|) arises from the regular sector, and the other 
one from the irregular sector, therefore 



p a {w) = Preg(u) + Pirr{u). 



(64) 



For comparison, we consider, first, the non-interacting 
limit, since in the absence of interactions simple ex- 
pressions can be obtained for the DOS. In this limit 
eg — 0, and the DOS associated with the ir- 
regular sector does not couple to the gapless fields. It is 
given by the DOS of the gapped system only, 



p^(w,x = 0) = ^-i?(M-A) T = 

2TTV y/u 



(65) 



Since p™ g (uj,x = 0) = 1/2ttvf, the total DOS exhibits a 
pseudo-gap for |w| < A corresponding to the pinning of 
half of the degrees of freedom. 

Let us now consider the general interacting case. In 
Sec. MI C II we have seen that the Green's function of 
the regular sector coincides with the Green's function 
of a spinlcss LL with the substitutions v — > vsll and 
7 — > Ireg- Hence, the corresponding DOS can be read off 
from Eq. (|61[) and is given by 



Preg(u) = 



l^/A-| 2 ^ 
2ttv S llF(1 + 2^ reg ) ' 



(66) 



where A~ = vsll/o, = V-/a. Note that the prefactor 2 
in Eq. (|61j) is absent because p reg refers here to a single 
mode (L I) or (R f). 

The irregular part of the DOS can, following Ref. Ill6l 
be expressed by 



delm[F+i ret (0, e)}lm[FZ{ ret (0, u - e) 



Pi rr (uj) 

-ret ^ 

which is proved in AppcndixICl The function F_± re (x, ui) 
is given by Eq. (fBU)) . The only part that is not purely 
real and contributes to Im[FTi re (0, oj — e)] is the part 
depending on the Bessel function Kx_ lirr . Expanding 

the Bessel function for small x (see Appendix IB"]) leads to 



A a T(2~/ irr ) 



Aa 



2 liT 



(68) 



The function F^ 1 (x, iu) n ) was given in Eq. (|55j) in form 
of a convolution integral and so is, in principle, more 
difficult to evaluate. However, for the DOS only the x — > 
limit is required, for which we can solve the Fourier 
integral for F_ t (x = 0,oj„) directly. Indeed, at x = 0, 
the Fourier transformation of Eq. (|5ip reduces to the 



standard integral of the Gamma function and we obtain 

i 

1 



'mv + 



r l- 



K 



fmv + \~ fmv + y 

iw„ + h 4W„ 

V 2K n ) \2K n J 



(69) 



2K ) V 2K 

With the analytic continuation iu n — > u+ we then find 
Im[F+i ret (0,w)] = -sign(w)tf {\uj\ - A) 

7T f\uj\-A 



A« + r(^; 



A+ 



(70) 



where we have used sin(^)F(l — -4=) = 7r/F(^) and 
defined A+ = v + /a as well as A = mv + /2K = 



B c ff At /AK. From the step function $ in the latter 
equation follows that pi rr (uj) = for \u>\ < A. For u> > A 
we then have 

1 

Pi rr (uj) 



A 



7ra (A-)27_(A+)^r(2 7irr )r(i) 

de(uj - e ) 2lirr ~ 1 (e - A)^ _1 , 

1 

7ra(A«)27^(A+)^r(2 7irr )r(i) 



(71) 



with B the Beta function. This result can be sim- 
plified by using 2^ irr + ± = 2K \ LL and B(x,y) = 
T(x)T(y)/T(x + y). Extending it furthermore to all lu 
we find 



0(|w|-A) (M-A) 31 ^' 



r( 



2 A -, 



na(Aa) 2 ^(A+)i< 



(72) 



Note that even though pi rr contains the gapless part 
ImFZi the threshold behavior of the massive sector oc 
$(|w| — A) is preserved in pi rr . The power-law at \uj\ > A, 
however, is the combination of the massive and the gap- 
less contributions and curiously turns out to be given by 
the LL parameter K$ll of the gapless fields (j>-,6- only, 
which is a behavior found also for the DOS of the gapped 
sector in charge density wave systems^ 

In the non-interacting limit, Ksll — > 1 and — > 0, 
the frequency dependence of p™. in the vicinity of \oj\ ~ 
A is recovered. However, there is a difference to Eq. (f6"5")) 
by a numerical prefactor of order 1. This is an artifact 
from the high-energy cutoff procedures and just reflects 
the unavoidable uncertainties with numerical prefactors 
in LL theory. 

The total DOS is subsequently obtained as the sum 
p(uj) = pf(u>) + pi(u) = 2[p reg (uj) +PirrM], where the 2 
arises from summing over the two regular {(L 4-), (R t)} 
as well as the two irregular {(L f), (R 4-)} contributions 
to the DOS. 
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FIG. 1: The local DOS as a function of energy for the SLL 
appearing through nuclear spin order in GaAs quantum wires 
(solid line). The LL parameters are chosen as2£ K c = 0.5 and 
K 3 — 1. Further system parameters as well as the size of the 
pseudogap A (after renormalization by electron interactions) 
are determined according to Ref. [H. Below A, the DOS is 
described only by preg, while at ui > A the DOS is dominated 
by the singular behavior of pi rr . For comparison, the dashed 
line shows the DOS of a regular spinful LL with the same 
parameters K c and K 3 . 



Let us note that since cross-correlators in the irregular 
sector between (L f) and (R I) modes are nonzero, the 
spin-polarization densities of state cr x (ui) or (J v {ui) probe 
directly pi rr , as they are linear combinations of (ijjt'4'l) 

and (ip\ipf). This is a consequence of the partial electron 
ordering in form of a spiral electron spin density wave, 
pinned to the nuclear spin spiral, or of the transverse 
ferromagnetic electron order aligned with the external 
magnetic field for the quantum wires with spin-orbit in- 
teraction at commensurability. 

In Figs. [1] and [5] we plot the DOS as a function of 
energy for the SLL appearing through nuclear spin or- 
dering in typical GaAs quantum wires^ for K s = 1 and 
K c > 1/3 as well as K c < 1/3. For K s = 1 the value 
K c = 1/3 marks the point at which Ksll = 1/2, and 
at which the exponent in Eq. (|72j) changes its sign. For 
K c > 1/3 the irregular contribution diverges at \u\ = A, 
while for K c < 1/3 it goes to zero. In both cases, how- 
ever, Figs. [T] and [2] show that pi rr dominates the DOS at 
\uj\ > A. A similar effect is observed when calculating the 
tunneling current between a LL and a superconductor 
which is not surprising, since both calculations involve 
convolutions between a LL-type power-law decaying com- 
ponent, and a gapped component. It needs to be stressed 
furthermore that our analysis is valid only in the vicin- 
ity of A, as further fluctuations of the gapped fields are 
not taken into account. However, as shown in a very 
recent complementary study of the SLL,— the inclusion 
of such fluctuations does not change the overall behavior 
of the DOS close to A but leads to modified power-law 
exponents. 



FIG. 2: Same as in Fig. [T] but for the reduced value K c = 0.3. 
Since K c < 1/3, the power-law singularity of pi rr at u> — A 
is replaced by pi rr — > 0. Yet, pi rr still dominates the DOS 
at uj > A. Note that A takes here a different value than in 
Fig. [T]because of the different renormalization by the stronger 
interactions.— 



V. CONCLUSIONS 

In this paper we have computed the position-frequency 
dependent Green's functions and the local DOS for vari- 
ous LLs. Most of the results are given in closed analytical 
form, to our knowledge, for the first time. These Green's 
functions are necessary for any local spectral characteri- 
zation of the LLs, and, in particular, we have used them 
to derive analytic expressions for the local DOS of all the 
investigated LLs. While the latter is a well-known result 
for the standard LLs ; 40 ' 41 we have put special attention 
to the HLL appearing at the edge of topological insula- 
tors and the SLL that appears in systems that are often 
used to mimic HLL behavior. We have shown that while 
HLL and SLL arc indeed equivalent for non-interacting 
electrons, they have substantially different spectral prop- 
erties for interacting electrons. While the HLL behaves 
as a standard spinless LL, the SLL shows with the ap- 
pearance of the irregular parts of the Green's functions 
entirely different response functions. The main reason 
for this striking difference is that the SLL exhibits the 
response of a standard gapless LL together with that of 
a gapped system with a dynamically created gap. The 
response therefore exhibits a pseudogap. Close to the 
pseudogap, the strong fluctuations of the fields conjugate 
to the gapped fields lead to the strong irregular response 
that dominates the spectral properties. Below the pseu- 
dogap, the local response described by the DOS shows a 
regular LL behavior, yet with modified exponents. Inter- 
estingly, the irregular behavior has an effect even below 
the pseudogap when considering the nonlocal response 
described by the Green's functions at i / 0. This be- 
havior is entirely absent in any HLL. Hence any analysis 
of the system properties and application that is based on 
the apparent equivalence of HLL and SLL needs to take 
this important difference into account. 

Concerning the detection of the discussed SLL prop- 
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erties it should be stressed that the SLL physics arises 
only when the chemical potential is kept near the center 
of the pseudogap. This occurs naturally in the nuclear- 
spin-ordered phase, in which the gap opens symmetri- 
cally about the chemical potential. It requires, however, 
a fine tuning of the electron density to kp = k so for the 
spin-orbit split quantum wires. This means, it is not pos- 
sible to probe SLL physics by modulating the chemical 
potential, and techniques such as STM need to be used. 

The systems that can exhibit the SLL physics are ID 
conductors in the nuclear ordered phase such as GaAs 
quantum wires^ or carbon nanotube o 37 ' 38 grown entirely 
from 13 C isotopes to provide the nuclear spins .^^r- 1 - 1 ^ In 
these systems, the particular spectral properties of the 
SLL could be used to indirectly prove the existence of 
the nuclear spin order. Furthermore, SLL physics is ex- 
pected to arise in ID conductors with strong spin-orbit 
interaction such as GaAs quantum wires, InAs nanowires, 
or Ge/Si core/shell nanowires in the presence of an exter- 
nal magnetic field ^ 13 ' 102 ' '-^ as well as carbon nanotubes 
in strong electric fieldsi 14 ' 15 



does not lie on a branch cut. From Eq. (|35j) we have 
-la 



2iza 



2ira 



x + v(t =F «0) 
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rx + v(t =p iO) 



x - v(t T *0)_ 

1 7+1 



It 



rx — v(t =F iO) 



(Al) 



with 7 = min(7i,7R) = {K+K~ l -2)/A. Since 7+1 > 1, 
the singularity at rx — vt = requires regularization, 
which can be done in the same way as in Eq. (|37|) by 
reducing the exponent to 7 by derivation with respect to 
-j^{d rx — d v t). We apply this directly on the retarded 
Green's function G r r e *{x,t) = d(t)[G> r (x, t) — G< r ], and 
obtain for the Fourier transform 



G r r f{x,oj) 



47T7 

with uj + = lu + iO and 



d rx + ^± 
v 



I(x,lu), (A2) 
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Appendix A: Retarded Green's functions in (x,co) 
space for the regular spinless Luttinger liquid 

In this appendix we provide the details of the calcu- 
lation of the real-time Fourier transform of the retarded 
Green's functions. 

The greater and lesser functions can be inferred di- 
rectly from the imaginary time Green's function. Indeed, 
since the latter is given by power- laws in x±ivr, the real 
time Green's functions can be obtained by a simple ana- 
lytic continuation of the time variable ; 41 ' 121 G(x, r — > 
i(t - iO)) = -iG>(x,t), and G(x,t -> i(t + iO)) = 
iG < (x,t). Due to the power-law forms, this substitu- 
tion is equivalent to the usual Wick rotation in frequency 
space and is valid as long as the analytical continuation 



I{x,u)) 



dte iuj+t 



x 2 +v 2 (t~i0y 



x 2 +v 2 (t + iO) 2 



(A3) 



The term depending on (t — iO) has branch cuts starting 
at ±|x|/u + iO running to ±00 slightly above the real 
axis, and the term depending on (t + iO) has the same 
branch cuts shifted slightly below the real axis [see Fig. 
[3] (a)]. Since to is real, we have I(x, — u>) = —I*(x,u>), 
allowing us to focus on u > 0. We can then deform the 
integration contour as shown in Fig. |3] (b). The only 
nonzero contribution comes from the two straight lines 
above and below the branch cut from \x\/v to +00. This 
leads to (with y = t + \x\/v) 

1 ' {v/a) 2 i J ' y [y(y + 2\x\/v)p 

(A4) 

Due to the -MO in w+ this integral is convergent and 
very similar to Eq. (|38|) . In fact, if the integration con- 
tour in Eq. (|38[) is deformed such that it runs around a 
branch cut, the integral becomes identical to Eq. (|A4[) 
under the replacement la; ^! — » — iuj+. This integral can 
be evaluated using Ref. Ill5l 3.383.8, and produces a 
result equivalent to Eq. ([3U1) . 



I(x, ui) 



2eu 



2i|x|w 



uT(7) \ us + a 



K i(\x\w + /iv). 



(A5) 

This form also holds for u) < as it satisfies I(x,u) = 
—I*(x,—oj). The branch as a function of lu runs along 
the negative imaginary axis. Therefore, the Wick rota- 
tion w+ — > iuj n is indeed problematic for uj n < 0, which 
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(a) 



table, Ref. I115L 3.387.7, we see that it evaluates to 



(b) 



\x\/v 



\x,uj) =2i dfe-H* 
Jo 



1 T 



v 2 t 2 



% — r(i-7) 



2xv 



Hi 



Yi 



(A7) 



with H Q the Struve function and Y a the Bcssel function 
of the second kind. This function, in combination with 
the solution I(x,u>) of Eq. (|A5[) . provides the result of 
Eq. ([SH]). 



FIG. 3: Branch cuts and integration contours for Eq. (IA3|) in 
the complex t plane, (a) The first term in Eq. (|A3f) has the 
upper right and lower left branch cuts, and the second term 
the upper left and lower right branch cuts. The integration 
runs along the real axis from t = to oo. (b) Since u > 0, 
the integrand converges in the upper half plane and can be 
deformed to a part running along the positive imaginary axis, 
and a part running around the upper right branch cut. Along 
the imaginary axis, both terms in the integrand are identical 
and cancel each other. The contribution from the circular 
contour around the branch point at \x\/v + iO vanishes. The 
two legs above and below the branch cut give Eq. (|A4|) . where 
— 2zsin(7T7) captures the phase difference between the upper 
and lower legs. 



explains why separate calculations were necessary for the 
real and imaginary time Green's functions. 

The final formula for the retarded Green's function is 
obtained in the same way as in Sec. IIII Al using Eq. (|40[). 
We find 



G r r e ^{x,tu) 



_ e irk F x fl / 2i\ x \v 



2 v ^Fr(7 + l) v 2 \w + b' 



K^_d\x\u)+/iv) - svga.(rx)K^ + i{\x\u)+/iv) 



Appendix B: Details for the calculation of the 
density of states 

We evaluate in this appendix the limit x — > for the re- 
tarded Green's function required for the DOS p{w). Let- 
ting x —¥ and rescaling Eq. (|A1[) or (|35|) by t —> cut 
we see that p{uj) oc oj 2g for the regular LL. However, we 
see also that taking the x — > limit before evaluating 
the Fourier integral over t makes the integral divergent. 
Taking the limit x — > after integration does not cure 
this problem, and the power-law divergences that arise 
from the extrapolation of the long- wavelength LL theory 
into the short distance x — > regime are still present. 
This is, however, an artifact of the lack of small distance 
cutoff in the model and should have no influence on the 
DOS; thus, these terms can be safely neglected. In fact, 
they can be entirely suppressed by the regularization via 
a further differentiation of the integrand as in Eq. (|37|) . 

The goal of the present appendix is to evaluate the 
prefactors of the w 27 terms for a typical LL. This analysis 
is based on an expansion of t he Bessel functions K a at 
small argument (see Ref. Illfjl . 9.6.9), 



K a {z) 



>o 1 
2 



r(a)(f)" + r(-a)(! 



{l + 0(z 2 )). 

(Bl) 

Introducing this expansion into Eq. (|42[) for a regular 
spinlcss LL and neglecting the divergent parts, we obtain 



(A6) 



40iT(l + 7) 



(B2) 



so 



which is the result reported in Eq. (|42l) . 

To finish this section, let us calculate the contribution 
of the vertical contour on the positive imaginary axis 
in Fig. (b). While this integral has dropped out for 
I(x, to), it actually contributes to the Green's function of 
the irregular sector of Eq. (|56|) . which is of the form of 
Eq. (|A3[) . yet with a + sign between the two power laws. 
We denote this contribution by I vort . Using the standard 



Pr(w) = ImG^(a:^0,w) 

1 cos(7T7)r(i — 7) 
~ 2^u 2 27 v^Fr(7 + 1) 



uia 



27 



(B3) 



Using the doubling theorem of Gamma functions [Ref. 

GH 6.1.18], T(2z) = (27r)- 1 / 2 2 2z -5r(z)r(z + i) with 
, _ 1 



7, we can further simplify this result to 



p r (uj) = 



1 



1 



2-kv r(i + 27) 



2l 



(B4) 
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FIG. 4: Integration contour C for Eq. (|C2[I . The black dots 
indicate the Matsubara poles. The wiggled lines mark the 
branch cuts at Im(z) = separating retarded and advanced 
parts of F+, and at Im(z) = ito n separating retarded and 
advanced parts of F-. 



Appendix C: Proof of Eq. (|67]) 

In this appendix we derive the identity given by Eq. 
(p7|) . To this end, it is useful to change to a finite temper- 
ature description with discrete Matsubara frequencies. In 
this case, the Green's function of Eq. (|T7)) has the generic 
form 



G(iw n ) = F + (icj m )F_ {iw n 



(CI) 



where we omit all prefactors and indices that are not 
necessary for the current calculation. We are here con- 
fronted with an ambiguity of how to choose the Matsub- 
ara frequencies in the bosonization formulation, as iui n 
should in principle be a fermionic frequency. Indepen- 
dently of the choice of ioj m one of the F± would then de- 
pend on a bosonic and the other one on a fermionic Mat- 
subara frequency In the limit of zero temperature, how- 
ever, all Matsubara frequencies become continuous, and 
so we make here the most convenient choice of bosonic 
iw n = 2nir I /3 and fermionic iuj m = (2m + l)w/ j3, for inte- 
gers n and m, and /? the inverse temperature. This choice 
guarantees an equivalent treatment of both F± and al- 
lows us to avoid unnecessary complications with poles on 
the branch cuts. In a standard way, we can replace the 
Matsubara sum by the integral 



G(iuj n ) 



c £ tanh (^ 



with C a contour surrounding all Matsubara frequencies 
iuj m . Since F + (z) has a branch cut at Im(z) = and 
F-(iu n — z) has a branch cut at Im(z) = iu> n , the con- 
tour C takes the form shown in Fig. |U Each branch cut 
separates the retarded and advanced parts of the corre- 
sponding F±. Since iuj n is bosonic, there is no pole on 
the branch cuts, and the contour integral reduces to the 
evaluation along the cuts, 



d{lLJn) = /I ^ tanh (-r) [F + t{uj/) F + dV{u},)] 



x F r _ et {iw n -w') + / — tanh. 

J-oo 27r * V 2 

x F r + et {iuj n + oj') [F* dv {-uj') - F r _ et {-Lo')} . (C3) 



We now perform the analytic continuation iuj n — > lu + = 
uj + 10. Here we need to emphasize that the analytic 
continuation is ill defined for the tanh, as it is periodic 
in iuj n . To maintain the uniqueness of the analytic con- 
tinuation, we must first drop the iui n dependence of the 
tanh and then perform the analytic continuation for the 
functions F± only. The imaginary part of the resulting 
retarded Green's function can then be written as (after 
shifting the integration variable in the second integral by 



ImG ret (w) 



tanh 



— tanh 



X Im[Fi e V)]Im[f: et (u; - w ')]. 



(C4) 



At zero temperature f3 — > oo and the difference of tanh 
restricts the integration to the interval (0,uj) for uj > 
or (cj,0) for to < 0, and provides a factor 2sign(w). The 
result is 



ImG ret (u) 



- / dLj'Im[F: et (u')]lm[FL et {o 
77 Jo 



/)], 
(C5) 



which corresponds to Eq. (JB7] 



1 M. Konig, H. Buhmann, L. W. Molenkamp, T. Hughes, 
C.-X. Liu, X.-L. Qi, and S.-C. Zhang, J. Phys. Soc. Japan, 
77, 031007 (2008). 

2 M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 



3 R. Egger, A. Zazunov, and A. Levy Yeyati, Phys. Rev. 
Lett. 105, 136403 (2010). 

4 D. Kong, J. C. Randel, H. Peng, J. J. Cha, S. Meister, K. 
Lai, Y. Chen, Z.-X. Shen, H. C. Manoharan, and Y. Cui 
Nano Lett. 10, 329 (2010). 



16 



5 H. Peng, K. Lai, D. Kong, S. Meister, Y. Chen, X.-L. Qi, 
S.-C. Zhang, Z.-X. Shen, and Y. Cui, Nature Mater. 9, 
225 (2010). 

6 P. Stfeda and P. Seba, Phys. Rev. Lett. 90, 256601 (2003). 

7 Y. V. Pershin, J. A. Nesteroff, and V. Privman, Phys. 
Rev. B 69, 121306(R) (2004). 

8 P. Devillard, A. Crepieux, K. I. Imura, and T. Martin, 
Phys. Rev. B 72, 041309(R) (2005). 

9 S. Zhang, R. Liang, E. Zhang, L. Zhang, and Y. Liu, Phys. 
Rev. B 73, 155316 (2006). 

10 D. Sanchez, L. Serra, and M.-S. Choi, Phys. Rev. B 77, 
035315 (2008). 

11 J. E. Birkholz and V. Meden, Phys. Rev. B 79, 085420 
(2009). 

12 C. H. L. Quay, T. L. Hughes, J. A. Sulpizio, L. N. Pfeiffer, 
K. W. Baldwin, K. W. West, D. Goldhaber-Gordon, and 
R. de Picciotto, Nature Phys. 6, 336 (2010). 

13 B. Braunecker, G. I. Japaridze, J. Klinovaja, and D. Loss, 
Phys. Rev. B 82, 045127 (2010). 

14 J. Klinovaja, M. J. Schmidt, B. Braunecker, and D. Loss, 
Phys. Rev. Lett. 106, 156809 (2011). 

15 J. Klinovaja, M. J. Schmidt, B. Braunecker, and D. Loss, 
Phys. Rev. B 84, 085452 (2011). 

16 L. Fu and C.L. Kane, Phys. Rev. Lett. 100, 096407 

(2008) . 

17 J. Alicea, Phys. Rev. B 81, 125318 (2010). 

18 J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, 
Phys. Rev. Lett. 104, 040502 (2010). 

19 J. D. Sau, S. Tewari, R. M. Lutchyn, T. D. Stanescu, and 
S. Das Sarma, Phys. Rev. B 82, 214509 (2010). 

20 S. Tewari, T. D. Stanescu, J. D. Sau, S. Das Sarma, New 
J. Phys. 13, 06500 (2011). 

21 R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. 
Lett. 105, 077001 (2010). 

22 Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 
105, 177002 (2010). 

23 J. Alicea, Y. Oreg, G. Refael, F. von Oppen, and M. P. 

A. Fisher, Nature Phys. 7, 412 (2011). 

24 J. Klinovaja, S. Gangadharaiah, and D. Loss, 
larXiv:1201. 01591 

25 AY. Kitaev, Phys. Usp. 44, 131 (2001). 

26 C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. 
Das Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

27 F. Hassler, A. R. Akhmerov, C.-Y. Hou, and C. W. J. 
Beenakker, New J. Phys. 12, 125002 (2010). 

28 M. Wimmer, A. R. Akhmerov, M. V. Medvedyeva, J. 
Tworzydlo, and C. W. J. Beenakker, Phys. Rev. Lett. 
105, 046803 (2010). 

29 M. Duckheim and P.W. Brouwer, Phys. Rev. B 83, 054513 
(2011). 

30 C. Y. Hou, E. A. Kim, and C. Chamon, Phys. Rev. Lett. 
102, 076602 (2009). 

31 A. Strom and H. Johannesson, Phys. Rev. Lett. 102, 
096806 (2009). 

32 J. C. Y. Teo and C. L. Kane, Phys. Rev. B 79, 235321 

(2009) . 

33 S. Das and S. Rao, Phys. Rev. Lett. 106, 236403 (2011). 

34 R. Egger and H. Schoeller, Phys. Rev. B 54, 16337 (1996). 

35 P. Simon and D. Loss, Phys. Rev. Lett. 98, 156401 (2007). 

36 P. Simon, B. Braunecker, and D. Loss Phys. Rev. B 77, 
045108 (2008). 

B. Braunecker, P. Simon, and D. Loss, Phys. Rev. Lett. 
102, 116403 (2009). 

38 B. Braunecker, P. Simon, and D. Loss, Phys. Rev. B 80, 



165119 (2009). 

39 D. Loss, F. L. Pedrocchi, and A. J. Leggett, Phys. Rev. 
Lett. 107, 107201 (2011). 

4(1 A. O. Gogolin, A. A. Nersesyan, and A. M. Tsve- 
lik, Bosonization and Strongly Correlated Systems, Cam- 
bridge Univ. Press, Cambridge, 1998. 

41 T. Giamarchi, Quantum Physics in One Dimension, Ox- 
ford Univ. Press, Oxford, 2005. 

42 O. M. Auslaender, A. Yacoby, R. de Picciotto, K. W. 
Baldwin, L. N. Pfeiffer, and K. W. West, Science 295, 
825 (2002). 

43 O. M. Auslaender, H. Steinberg, A. Yacoby, Y. 
Tserkovnyak, B. I. Halperin, K. W. Baldwin, L. N. Pfeif- 
fer, and K. W. West, Science 308, 88 (2005). 

44 H. Steinberg, G. Barak, A. Yacoby, L. N. Pheiffer, K. W. 
West, B. I. Halperin, and K. L. Hur, Nature Phys. 4, 116 
(2008). 

45 Y. Jompol, C. J. B. Ford, J. P. Griffiths, I. Farrer, G. A. 
C. Jones, D. Anderson, D. A. Ritchie, T. W. Silk, and A. 
J. Schofield, Science 325, 597 (2009). 

46 M. Bockrath, D. H. Cobden, J. Lu, A. G. Rinzler, R. E. 
Smalley, L. Balents, and P. L. McEuen, Nature 397, 598 
(1999). 

47 Z. Yao, H. Postma, L. Balents, and C. Dekker, Nature 
402, 273 (1999). 

48 A. Bachtold, M. de Jonge, K. Grove-Rasmussen, and P. 
L. McEuen, M. Buitelaar, and C. Schonenberger, Phys. 
Rev. Lett. 87, 166801 (2001). 

49 Y. Ihara, P. Wzietek, H. Alloul, M. H. Riimmeli, Th. 
Pichler, and F. Simon, EPL 90, 17004 (2010). 

50 A. M. Chang, L. N. Pfeiffer, and K. W. West, Phys. Rev. 
Lett. 77, 2538 (1996). 

51 E. Slot, M. A. Hoist, H. S. J. van der Zant, and S. V. 
Zaitsev-Zotov, Phys. Rev. Lett. 93, 176602 (2004). 

52 A. N. Aleshin, H. J. Lee, Y. W. Park, and K. Akagi, Phys. 
Rev. Lett. 93, 196601 (2004). 

53 L. Venkataraman, Y. S. Hong, and P. Kim, Phys. Rev. 
Lett. 96, 076601 (2006). 

54 J. D. Yuen, R. Menon, N. E. Coates, E. B. Namdas, S. 
Cho, S. T. Hannahs, D. Moses, and A. J. Heeger, Nature 
Mater. 8, 572 (2009). 

55 C. Blumenstein, J. Schaefer, S. Mietke, A. Dollinger, M. 
Lochner, X. Y. Cui, L. Patthey, R. Matzdorf, and R. 
Claessen, Nature Phys. 7, 776 (2011). 

56 A. V. Moroz, K. V. Samokhin, and C. H. W. Barnes, 
Phys. Rev. Lett. 84, 4164 (2000). 

57 A. V. Moroz, K. V. Samokhin, and C. H. W. Barnes, 
Phys. Rev. B 62, 16900 (2000). 

58 M. Governale and U. Ziilicke, Phys. Rev. B 66, 073311 
(2002). 

59 A. Iucci, Phys. Rev. B 68, 075107 (2003). 

60 Y. Yu, Y. Wen, J. Li, Z. Su, and S. T. Chui, Phys. Rev. 
B 69, 153307 (2004). 

61 V. Gritsev, G. I. Japaridze, M. Pletyukhov, and D. 
Baeriswyl, Phys. Rev. Lett. 94, 137207 (2005). 

62 F. Cheng and G. Zhou, J. Phys.: Condens. Matter 19, 
136215 (2007). 

63 J. Sun, S. Gangadharaiah, and O. A. Starykh, Phys. Rev. 
Lett. 98, 126408 (2007). 

64 S. Gangadharaiah, J. Sun, and O. A. Starykh, Phys. Rev. 
B 78, 054436 (2008). 

65 A. Schulz, A. De Martino, P. Ingenhoven, and R. Egger, 
Phys. Rev. B 79, 205432 (2009). 

66 G. I. Japaridze, H. Johannesson, and A. Ferraz, Phys. 



17 



Rev. B 80, 041308(R) (2009). 

67 A. Schulz, A. De Martino, and R. Egger, Phys. Rev. B 
82, 033407 (2010). 

68 A. Strom, H. Johannesson, and G. I. Japaridze, Phys. 
Rev. Lett. 104, 256804 (2010). 

69 M. Malard, I. Grusha, G. I. Japaridze, and H. Johannes- 
son, Phys. Rev. B 84, 075466 (2011). 

70 C. Wu, B. A. Bernevig, and S.-C. Zhang, Phys. Rev. Lett. 
96, 106401 (2006). 

71 S. Gangadharaiah, B. Braunecker, P. Simon, and D. Loss, 
Phys. Rev. Lett. 107, 036801 (2011). 

72 E. M. Stoudenmire, J. Alicea, O. A. Starykh, and M. P. 
A. Fisher, Phys. Rev. B 84, 014503 (2011). 

73 E. Sela, A. Altland, and A. Rosch, Phys. Rev. B 84, 
085114 (2011). 

74 R. M. Lutchyn and M. P. A. Fisher, Phys. Rev. B 84, 
214528 (2011). 

75 S. Eggert and I. Affleck, Phys. Rev. B 46, 10866 (1992). 

76 E. Wang and I. Affleck, Nucl. Phys. B 417, 403 (1994). 

77 S. Eggert and I. Affleck, Phys. Rev. Lett. 75, 934 (1995). 

78 M. Fabrizio and A. O. Gogolin, Phys. Rev. B 51 17827 
(1995). 

79 S. Eggert, H. Johannesson, and A. E. Mattsson, Phys. 
Rev. Lett. 76, 1505 (1996). 

80 S. Eggert, A. E. Mattsson, and J. M. Kinaret, Phys. Rev. 
B 56, R15537 (1997). 

81 A. E. Mattsson, S. Eggert, and H. Johannesson, Phys. 
Rev. B 56, 15615 (1997). 

82 V. Meden and K. Schonhammer, Phys. Rev. B 46, 15753 
(1992). 

83 J. Voit, Phys. Rev. B 47, 6740 (1993). 

84 S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, 
J. M. Tranquada, A. Kapitulnik, and C. Howald, Rev. 
Mod. Phys. 75, 1201 (2003). 

85 P. Kakashvili, H. Johannesson, S. Eggert, Phys. Rev. B 
74, 085114 (2006). 

86 I. Schneider and S. Eggert, Phys. Rev. Lett. 104, 036402 
(2010). 

87 D. Schuricht, F. H. L. Essler, A. Jaefari, and E. Fradkin, 
Phys. Rev. Lett. 101, 086403 (2008). 

88 D. Schuricht, F. H. L. Essler, A. Jaefari, and E. Fradkin, 
Phys. Rev. B 83, 035111 (2011). 

89 D. Schuricht, Phys. Rev. B 84 , 045122 (2011). 

90 D. Schuricht. [arXiv:1112. 30451 

91 M. I. D'yakonov and A. V. Khaetsku, JETP Lett. 33, 110 
(1981). 

92 M. I. D'yakonov and A. V. Khaetskii, Sov. Phys. JETP 
55, 917 (1982). 

93 O. A. Pankratov, S. V. Pakhomov, and B. A. Volkov, 
Solid State. Comm. 61, 93 (1987). 

94 M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buh- 
mann, L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Sci- 

^ ence 318, 766 (2007). 

95 A. Roth, C. Briine, H. Buhmann, L. W. Molenkamp, J. 
Maciejko, X.-L. Qi, and S.-C. Zhang, Science 325, 294 
(2009). 

96 C. Xu and J. E. Moore, Phys. Rev. B 73, 045322 (2006). 

97 Y. Tanaka and N. Nagaosa, Phys. Rev. Lett. 103, 166403 
(2009). 



us 

119 



C.-X. Liu, J. C. Budich, P. Recher, and B. Trauzettel, 
Phys. Rev. B 83, 035407 (2011). 

M. Hohenadler, T. C. Lang, and F. F. Assaad, Phys. Rev. 

Lett. 106, 100403 (2011). 

M. Hohenadler and F. F. Assaad. [arXiv:1110.3322l 
Note that the signs in Eqs. (|15[) and (|16|) correct the typo 
in the signs of Eqs. (45) and (46) of Ref. [3^ and agree 
with the supplementary online EPAPS document No. E- 
PRLTAO-102-025914 of Ref. H3- 

C. Kloeffel, M. Trif, and D. Loss, Phys. Rev. B 84, 195314 
(2011). 

X.-J. Hao, T. Tu, G. Cao, C. Zhou, H.-O. Li, G.-C. Guo, 
W. Y. Fung, Z. Ji, G.-P. Guo, and W. Lu, Nano Letters 
10, 2956 (2010). 

W. Lu, J. Xiang, B. P. Timko, Y. Wu, and C. M. Lieber, 
PNAS 102, 10046 (2005). 

J. Xiang, W. Lu, Y. Hu, Y. Wu, H. Yan, and C. M. Lieber, 
Nature 441, 489 (2006). 

J.-S. Park, B. Ryu, C.-Y. Moon, and K. J. Chang, Nano 
Letters 10, 116 (2010). 

H. Lee and H. J. Choi, Nano Letters 10, 2207 (2010). 
L. N. Pfeiffer, A. Yacoby, H. L. Stormer, K. L. Baldwin, 
J. Hasen, A. Pinczuk, W. Wegscheider, and K. W. West, 
Microelectron. J. 28, 817 (1997). 

F. Simon, C. Kramberger, R. Pfeiffer, H. Kuzmany, V. 
Zolyomi, J. Kurti, P. M. Singer, and H. Alloul, Phys. 
Rev. Lett. 95, 017401 (2005). 

M. H. Riimmeli, M. Loffler, C. Kramberger, F. Simon, F. 
Fiilop, O. Jost, R. Schonfelder, A. Griineis, T. Gemming, 
W. Pompe, B. Biichner, and T. Pichler, J. Phys. Chem. 
C 111, 4094 (2007). 

H. O. H. Churchill, A. J. Bestwick, J. W. Harlow, F. 
Kuemmeth, D. Marcos, C. H. Stwertka, S. K. Watson, 
and C. M. Marcus, Nature Phys. 5, 321 (2009). 

H. O. H. Churchill, F. Kuemmeth, J. W. Harlow, A. J. 
Bestwick, E. I. Rashba, K. Flensberg, C. H. Stwertka, T. 
Taychatanapat, S. K. Watson, and C. M. Marcus, Phys. 
Rev. Lett. 102, 166802 (2009). 

R. Egger and A. O. Gogolin, Phys. Rev. Lett. 79, 5082 
(1997). 

C. Kane, L. Balents, and M. P. A. Fisher, Phys. Rev. Lett. 
79, 5086 (1997). 

I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, 
Series, and Products, 5th edition, Academic Press, 1994. 
O. A. Starykh, D. L. Maslov, W. Hausler, and 
L. I. Glazman, in Low- dimensional Systems: Inter- 
actions and Transport Properties, Lecture Notes in 
Physics 544, T. Brandes (ed.), 37-78, Springer, 2000 
( |arXiv:cond-mat/9911286j ). 

J. Voit, Eur. Phys. J. B 5, 505 (1998). 
P. B. Wiegmann, Phys. Rev. B 59, 15705 (1999). 
M. Abramowitz and I. A. Stegun, Handbook of Mathe- 
matical Functions, Applied Mathematics Series, Vol. 55, 
National Bureau of Standards (Washington), tenth print- 
ing, 1972. 

C. Bena, Phys. Rev. B 82, 035312 (2010). 
J. W. Negele and H. Orland, Quantum Many-Particle Sys- 
tems, Frontiers in Physics, Addison- Wesley, Reading MA, 
1988. 



